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We describe how one can use chiral perturbation theory to obtain results for physical 
quantities, such as quark masses, using partially quenched simulations. 

For some time to come, simulations of lattice QCD will not work directly with 
physical up and down quarks. This is because the computer time required scales 
roughly as m"^'^ with existing algorithms (at fixed physical box size). Present 
simulations, using machines sustaining up to 0.3 Teraflops, are limited to quark 
masses greater than about half the strange quark mass, an order of magnitude larger 
than the average up and down quark mass. Even a dedicated machine sustaining 
10 Teraflops (hopefully to be available in 2003-4) will allow quark masses to be 
reduced to only ~ ms/8. 

Thus an extrapolation in light quark masses is required. Fortunately, this can 
be done, using chiral perturbation theory (ChPT), once the simulated masses are 
small enough. In practice, it is feasible to determine the functional forms with 
which to extrapolate at next-to- leading order (NLO) in the chiral expansion. A 
useful way of thinking about the extrapolation is that simulations with moderately 
light quark masses can be used to determine the parameters of the chiral Lagrangian 
(including the Gasser-Leutwyler coefficients Li-io which appear at NLO), and then 
the extrapolation can be done "by hand" . 

In this talk we describe some recent work in which we show how the c 
extrapolation can be aided by the use of partially quenched (PQ) simulationsEI 
These are simulations in which the "valence" and "sea" (= "dynamical") quarks 
have different masses. The key theoretical observation is that, if both valence and 
sea quarks are light enough, then the chiral Lagrangian describing the long distance 
properties of the PQ simulations contains the same parameters (/, {qq)^ and 
as appear in the chiral Lagrangian for QCD □. This follows from the work of Ref. 
Thus one can extrapolate to the physical theory using PQ as well as unquenched 
simulations, without introducing new, unphysical, parameters. This is true despite 
the fact that the PQ theory itself is unphysical. ° 

The practical importance of this observation has yet to be seen, but could be 
significant. It is relatively cheap (CPU ~ Triy^) to reduce the valence quark masses 



*Based on talks given by S. Sharpe and N. Shoresh at DPF 2000, August 2000. 

"A clear example of the unphysical nature of PQ theories is that some flavor-singlet correlation 

functions have double poles. 
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at fixed sea quark mass, and, in fact, most calculations including sea quarks do 
undertake extensive PQ simulations. Our point (amplified in examples given below) 
is that these simulations should not be viewed as a somewhat improved quenched 
approximation, from which one can only obtain qualitative information, but rather 
as a quantitative technical tool that one can use to obtain physical parameters. 

Note that it is essential for the PQ simulations to be done with the same number 
of light dynamical quarks as in physical QCD, since the parameters of the chiral La- 
grangian depend on this number. This brings up the tricky question of whether the 
strange quark is light enough to be described accurately enough by ChPT at NLO. 
The answer depends on the physical quantities being considered. For quantities in- 
volving the light pseudoscalar mesons, ChPT including the strange quark seems to 
be reasonably convergent [NLO corrections are generically ~ M|-/(47r/7r)^]. Thus it 
is probably appropriate to work in the sector with three light quarks. For baryonic 
quantities, where the expansion parameter is Mk / {^t^ f-w), the situation is much 
less clear. Here one may be forced to consider a two-flavor chiral Lagrangian with 
ms-dependent parameters, in order to make use of PQ simulations. In any case, the 
beauty of lattice simulations is that they can be used, in principle, to resolve these 
questions by comparing numerical results to the different theoretical predictions. 

We now turn to a more detailed explication of these general remarks. This will 
be done using the properties of pseudo-Goldstone mesons (henceforth referred to 
generically as "pions") as an example. More details can be found in Refs^Hu. 

We begin by recalling the form of the chiral Lagrangian for PQQCD B, 

£ = ^str {d^Ud^U^) - ^str {xU^ + Ux) - L, [str (xC/t + C/x)] ' + • ■ • ,(1) 

where, for brevity we have displayed only one of the NLO terms. Note that C 
has the same form as that for QCD, except that U = exp(2in//) is an element of 
the graded group SU{5\2) rather than SU{3) (assuming three sea quarks and two 
valence quarks), that traces are replaced by supertraces (str), and that the mass 
matrix which appears in x = 2^M (with /i = — (qq)/f'^) is enlarged to 

M= dia,g{mA,mB,mu,md,ms,mA,mB) (2) 

with mA,B being the masses of the valence quarks and the corresponding ghosts.** 
An important theoretical point is that the "pion" field n is "straceless" , strH = 
0, so that U is an element of SU{5\2) rather than C/(5|2). Since the symmetry 
breaking pattern is SU{5\2)l x SU{5\2)ii SU{5\2)v, this means that there is one 
element of H for each Goldstone particle. This all sounds very reasonable and hardly 
worthy of note. We mention it for two reasons. First, because it is in sharp contrast 
to the situation in the fully quenched theory. The quenched chiral Lagrangian must 
be constructed using a field which is an element of the unitary graded group, e.g. 
S G U{2\2) for two valence quarks. The non-anomalous chiral symmetries then 

''In the subsequent expressions, these masses can be taken itp be lattice quark masses, defined 
through Ward Identities, up to corrections of 0(a). See Ref. |3 for more discussion. 
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allow arbitrary functions of $o « InsdctE in the Lagrangian (as long as they are 
consistent with parity), including the well-known mass term ttiq^q. The theory is 
then no longer predictive, because there are an infinite number of new couplings, 
and no power counting scheme with which to order them. One must proceed in a 
phenomenological manner, assuming that most couplings are small, and ignoring 
higher-order $o loops. The essence of the problem is that the physics due to the 
scale Too ~ 1 GeV cannot be integrated out. 

The second reason we raise this point is that, if you look at papers doing 
calculations in PQChPT, including our own, you will find that U is enlarged to 
S = J7exp(i2$o/\/3/) € U{5\2). Given the discussion in the previous paragraph, 
you might well wonder why one would include the $0 lose predictive power. 
In our case, it was because we did not realize until recently that one could work 
without the $0, and so we followed a similar path to that for the quenched theory. 
The good news, however, is that one can show that if one takes previous results 
obtained including the $0i and sends TOo^— > 00, one recovers the results that would 
have been obtained if $0 = all along H. In fact, it is technically easier to do the 
calculations this way, rather than to explicitly project against $o- 

In summary, for light enough sea quark masses one can non-perturbativelj;, in- 
tegrate out the $0 from the PQ theory. This result was anticipated in Refs. Uu on 
the basis of a one-loop analysis. 

As our first application, we have calculated the pion masses and decay constants 
for PQ theories with an arbitrary number of sea quarks with arbitrary masses §, 
extending the work of Refs. 0i for degenerate sea quarks. This allows a determina- 
tion of /, fj, and Lis- We find, for example, that the mass of the fiavor non-singlet 
meson created by the operator Aj^B is 



Mis = Xv 



8N 8N 
1 + —{2U-Li)xs + -j^{2Ls-L5)xv + chiral logs + 0{p^) 



(3) 

where Xv s are the average valence and sea quark masses, respectively, multiplied 
by 2/i. The "chiral logs" are generically of the form c(x//^)lnx, with c a known 
number, and so do not introduce any new constants. The dependence on the scale in 
the logarithm is cancelled by the implicit scale dependence of the Gasser-Leutwyler 
coefficients. Obtaining these logarithms for arbitrary quark masses is the hard part 
of the calculation. Note that Eq. (||) includes results for the unquenched theory; 
these can be obtained by setting niA and TOb equal to sea quark masses. 

Fitting Eq. to numerical results (assuming / is known from fits to the decay 
constants) one can, in principle, determine 2LQ—L4 and 2L^—Lc,. With unquenched 
simulations alone, however, a separate determination of these two combinations re- 
quires non-degenerate sea quarks (for otherwise Xv = Xs)- contrast, with PQ 
simulations a separate determination is straightforward since Xs and Xv are inde- 
pendent, and one need only use degenerate sea quarks. It is particular noteworthy 
that the combination 2Lg,—L^, which is related to the value of the physical up-quark 
mass as discussed below, can be determined from the dependence on Xv^ using only 
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a single mass for the degenerate sea quarks Q. Of course, this mass must be hght 
enough that NNLO corrections are smaU. 

This discussion shows how the extra "dials" that one can adjust in the PQ theory 
simphfy the extrapolation to unquenched QCD. We stress, however, that we are not 
particularly advocating the use of three degenerate sea quarks — the so-called "2 + 1" 
simulations with rUu = ma < rus work just as well, and might be preferred as they 
approach closer to the physical parameters. Our general formulae for the chiral logs 
apply to either case. 

It seems that it will be practical to carry out such determinations in the next 
few years, and we spend the rest of this talk discussing issues that might arise. 
First recall that we only have good experimental information on two combinations 
of L4_8 (all values are quoted at the scale M^): 

L5 = (2.3±0.2) X 10"^ 2L7 + L8 = (0.4±0.1) X 10"^ (4) 

The remaining 3 linear combinations will either be very difficult, or impossible i, 
to learn from experiment alone. One can estimate these 3 combinations using 
phenomenological models, and there are a fairly standard set of values that can be 
found in texts and review articles. Thus, determining L4_g will allow both a test 
of QCD (Do these physical parameters equal the experimental values of Eq. (^)?) 
and a test of phenomenological models. 

Perhaps the most interesting aspect of any such determination is the light it 
might shed on the value of to„, the physical up-quark mass.'^ A vanishing value 
solves the strong-CP problem. At NLO in ChPT, and taking the experimental 
meson masses as inputs, the value of m,i depends on the combination iL^—L^. The 
standard parameters (which correspond to mulmd ~ 1/2) give 2Lg, — Lc, w 0, while 
to obtain m„ — one needs a value in the range (—1.2) — (—2.6) x 10"'^. There is a 
large uncertainty since, in this scenario, NLO corrections are « 50%, and so NNLO 
uncertainties are large. The upper and lower ends of this range come from writing 
the NLO expressions either in terms of quark masses or of squared meson masses. 
Because of this uncertainty, one can only hope to rule out m„ = by determining 
the Li. To demonstrate that m„ = will be much more difhcult. 

We have examined how results from PQ simulations would look in these two 
scenarios: the "Standard L's" or "m,i = L's" ."^ For simplicity, we set m„ = m^, and 
consider only charged meson properties, so L7 does not enter. We set L4 = Lq = 
in both scenarios, and^ 

{L5,Ls) X 10^ « (2.3, 1.2) [standard] or (1.2, 0) [to„ = 0] (5) 

Our values of / and the quark masses are then chosen to reproduce the physical 
ratios M^^/Mk, fTr/fK, and /jr/Af^, within a small tolerance. We stress that our 

'^This point has been stressed in Ref. S, where more details of the following discussion can be 
found. PI 

''Some results for the former appear in Ref. p, while those for the latter are new. 

'^For the m„ = scenario, L5 is about half the value quoted in Eq. (kh because we write the NLO 

expressions in terms of quark masses. 
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Standard L 's "ra„ = 0" L 's 
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Fig. 1. Contour plots of 5^. The x- and y-axes are the light sea and valence quark masses, 
respectively, in units of the physical strange quark mass. In the left-hand plot, the white region 
has < (5^ < 0.1, and the contours increase from 0.1 to 0.2 (upper left) in steps of 0.02. In the 
right-hand plot the contours start at 0.475 (lower left) and increase to 0.825 in steps of 0.025. 



choices of parameters are not unique — indeed, there is not enough experimental 
information to fix all the parameters — but rather that they are representative of 
the two scenarios. 

Meson properties depend on four parameters: the two valence quark masses and 
the light and strange sea quark masses. To reduce this to two parameters, we show 
results for "kaon-like" mesons, those in which the masses one of the valence quarks 
and the strange sea quark are fixed to the physical strange quark mass (mst). The 
properties then depend on the light valence and sea quark masses. If both are equal 
to the physical average light quark mass then our kaon is very close to the physical 
kaon. The quantities we plot are the ratios of NLO to LO contributions: 

4' = (A^/Xv) - 1 • and 4 = (/k//)-1. (6) 

These should have magnitudes small compared to unity for the chiral expansion to 
be reliable. The kaon-like mesons include some examples of the poorest convergence 
of the chiral expansion. 

In Fig. we show contour plots of 5^ for the two sets of parameters. Unquenched 
simulations lie on the diagonal, while most PQ simulations lie below the diagonal. 
Most present simulations lie to right of the x = 0.5 line. We observe that the chiral 
expansion is well behaved for standard parameters (and, indeed, the PQ region 
below the diagonal has better convergence than the unquenched line), while the 
expansion is not reliable for the m„ = case. All one can say in the latter case is 
that there should be substantially more curvature in the plots of versus quark 
mass than for standard parameters. 

In Fig. § we plot the NLO/LO ratio for decay constants of kaon-like mesons. 
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Standard L's 
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Fig. 2. Sj, along rays of fixed angle in the contour plots. 



We plot slices through the contour plot of Sj^ (not itself shown) along the various 
"rays" shown by the dashed lines in Fig. |l|, with the angle being that w.r.t. the 
X-axis. In this case the chiral expansion is reasonably well behaved for both sets of 
parameters, although the results themselves differ substantially. 

These plots bring out a warning to those fitting PQ or unquenched data. If 
one were to fit the curves for, say, 15° — 45°, in the region 0.5 < a; < 1 with a 
straight line and extrapolate, one would make a 10 — 20% error in the final result 
for fx- This is exactly the fit which would result if one kept only the analytic 
terms, and dropped the chiral logs, in expressions such as Eq. (||). Clearly one must 
keep the logarithms when extrapolating. We stress that this does not not introduce 
new parameters, and that this is true not only for PQ extrapolations but also for 
extrapolations along the unquenched line at 45°. 
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